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Q_| Phenonienological Tsallis fits to the CMS and ATLAS transverse spec- 

I , tra of charged particles were found to extend for pT from 0.5 to 181 GeV 

in pp collisions at LHC at ^/s =7 TeV, and for px from 0.5 to 31 GeV at 
(— I I ^/s =0.9 TeV. The simplicity of the Tsallis parametrization and the large 

range of the fitting transverse momentum raise questions on the physical 
meaning of the degrees of freedom that enter into the Tsallis distribution 
or g-statistics. 
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Recently there is a lot of interest in the Tsallis fit to the transverse 
momentum data of charged particles measured at very high energies of RHIC 
and LHC experiments [H O [3l HI [5l |6] . By this one understands that the 
use of the Tsallis distribution with a normalization constant Cq 

1 

hq{pT)=Cq l-(l-g)^ "\ (1) 



describes the experimental transverse momentum distribution data. The 
Tsallis distribution can be regarded as a generalization of the usual ex- 
ponential (Boltzmann-Gibbs) distribution, and converges to it when the 
parameter q tends to unity [7J, 

hq{pT)'^C^e^p(-P^y (2) 
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This approach is a nonextensive generahzation of the usual statistical me- 
chanics (characterized by a new parameter q, the nonextensivity parameter) 
and has been very successful in describing very different physical systems in 
terms of statistical approach, including multiparticle production processes 
at lower energie^. 

On the other hand, long time ago Hagedorn proposed the QCD in- 
spired empirical formula describing the data of the invariant cross section 
of hadrons as a function of pT over a wide range 




for pt — )• 0, 
for Pt — )• oo. 



(3) 



where C, po, and n are fitting parameters. This becomes a purely expo- 
nential function for small px and a purely power law function for large pT 
value^. It coincides with Eq. ([1]) for 

1 T 

n = and po = -. (4) 

q — l Q — 1 

Usually, both formulas are treated as equivalent from the point of view of 
phenomenological fits and are often used interchangeably [H El [3l [H El [6] . 
We follow this attitude for a while and shall discuss the possible physical 
implications later. 

For phenomenological as well as theoretical interests, it is useful to ex- 
plore where the Tsallis fit begins to fail at higher and higher px- Here we 
concentrate only on the recent high-py data of CMS [5j and ATLAS [3j Col- 
laborations at LHC. Excellent fit to the pT spectra was obtained there with 
the Tsallis and/or Hagedorn distributions for pj- from 0.5 GeV up to 6 GeV, 
in pp collisions at ^/s =7 TeV. However, CMS data extend to much higher 
range of pT, up to ~200 GeV/c. It will be interesting to know whether the 
good Tsallis fit continues to higher transverse momenta and how it would 
relate to results of fits to data at the lower energy of ^/s =0.9 TeV. 

The CMS [5J and ATLAS [Sj {Ed? N^^/ dp^)^ data shown in Fig. □ are 
taken essentially at the same kinematical windows, namely they correspond 
to an average over the data from rj=—r]Q to rj=+7]Q with r/o = 2.4 for the 
CMS measurements, and 2.5 for the ATLAS measurement^. To get the 



For a summary of earlier attempts of using Tsallis fits and detailed explanations of 
the possible meaning of the q parameter, together with up-to-date literature on this 
subject, see piQlfTO], 

Actually this QCD inspired formula was proposed earlier in [121 113) . 
We do not include here ALICE data [B] because they are for a smaller window, 
—0.8 < 77 < 0.8, and the data points are slightly higher for large pt because of the 
slight r] dependence of the spectra[B]. 



Tsallisfits printed on December 6, 2012 



3 



theoretical results from the distribution for comparison with experimental 
data we, therefore, calculate numerically 



E- 



ch 



1 r.4(E^' 

2% J-vo drj \ dp-^ 



(5) 



Here 



dy 
dr] 



rrij, cosh y 



(6) 



which is from Eq. (2.31) of [14j, and y is a function of rj and pT (Eq. 2.29) 
of 



■y = - In 
^ 2 



Pji cosh 1] + + pt sinh r] 



cosh rj + m? — px sinh 7] 



(7) 



To provide a theoretical fit to the experimental data, we follow the CMS 
Collaboration [1] and consider the differential cross section with a transverse 
Tsallis distribution [71 [15] in the form 



dp^ 



dy 



1 + ^ 

nT 



where 



Ej' 



m- 



+ pj. -m, 



(8) 



(9) 



and we assume m = m-,^ = 0.14 GeV. If we assume now a rapidity plateau 
structure with a constant CdN^h/dy, then the integral is 



E 



ch 



C dK 



ch 



dp-^ 



2r]o dy 



vo dy 
drj— 
-vo dr] 



1 + 



Ex 
nT 



(10) 



For each value oipx, the curves in Fig. [T]are obtained from such a numerical 
integration over rj. The px spectrum is therefore described by an overall 
constant A = CdN^h/dy and the parameters n and T. 

Previously, excellent Tsallis fit to the CMS px spectra at y/s = 7 TeV 
was obtained for 0.5 < px < Q GeV/c by using n = 6.6 and T = 0.145 
GeV [Ij. By using this set of parameters as initial guess, we search for the 
fits to the spectrum from 0.5 GeV up to 181 GeV. We find that for the set 
including the data points at the higher px region, the best fit is obtained 
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Tsallis fit for pp collisions 
q= (n + I) / II 
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Thermal disttibution 
(Tsallis disttibution with q-l) 
q= (n + 1)/ n 
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Fig. 1. (Color online) Left panel: Tsallis fits ^ to the CMS jSj and ATLAS 
[3] Collaborations data for pp at 7 and 0.9 TeV [5]. Right panel: the same data 
compared with the corresponding q — 1 (or n — >■ oo) curves. 



with A =4.06, n=6.6, and T = 0.147 GeV, which is essentially the same as 
the set obtained previously in [4J (Fig. [H left panel). What is surprising 
is that the Tsallis distribution that describes the data well at lower pT can 
describe the CMS data at high pT just as well. There is no departure of 
data from the Tsallis distribution from pT=0.5 GeV up to pT~200 GeV. For 
additional comparison, we include the ATLAS data in Fig. [TJ The ATLAS 
measurement has a slightly larger r] window, |t/|atlas ^ 2.5, instead of 
CMS's |r?|cMS < 2.4, and its pT values extends from 0.5 GeV/c to 36 GeV/c 
for ^=7 TeV. The ATLAS data at ^/s=7 TeV are consistent with the CMS 
data and the Tsallis fit (Fig. [H left panel). 

On the other hand, at the lower energy of ^/s = 0.9 TeV, we find that the 
pp data can be described well by the higher value of n=7.65 (with A=4.01 
and r=0.128 GeV) in FigH] Again, the CMS data and the ATLAS data are 
consistent with each other. The Tsallis distribution gives a good description 
of the V^=0.9 TeV data from pT=0.5 GeV/c to pt=31 GeV/c (Fig. [fl left 
panel) . 

It is instructive to visualize the difference between the Tsallis distribu- 
tion and the corresponding thermal distribution characterized by the same 
temperature T and by q = 1 (or n — )■ oo). This is shown on the right panel 
of Fig. [TJ Changing only the temperature parameter T would not give a 
better fit to the experimental data, to do this it is necessary to allow q to 
vary and become q > 1. 

The overall good agreement of the Tsallis parametrization and the ex- 
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perimental data is amazingly good. The absence of a departure of data from 
the Tsalhs distribution from px'=0.5 GeV to 181 GeV indicates that there 
are essentiahy only three degrees of freedom that count for the description of 
the pt distribution: an overall magnitude, and two other degrees of freedom 
to describe the shape. This result should be confronted with the (apparently 
equally successful) many parameter fits also presented in [Sj HI O |6] and us- 
ing known Monte Carlo programs. It can be interpreted as indication that, 
whereas hadronizing system formed in the process of particle production is 
very complex, only few degrees of freedom are really important. At lower 
energies (and for lower values of px) this was regarded as indication that 
such system can be described by simple (with q = 1) ox generalized (with 
q > 1) statistical models i8J. The results presented here indicate that, 
either such models work also for such large px or there are some specific, 
dynamical rather than purely statistical, phenomena at work as advocated 
recently in |16] . In any case, we need to take the physics contents of the 
Tsallis model seriously. 

At this point let us come back to formula (or pil|) ). It was proposed 
in [12] long time before Tsallis works [7j (and followed later by [13] and 
[11] ) with a simple aim: to phenomenologically interpolate between the soft 
region of pt— ^0, characterized by exponential behavior of pT distributions, 
and the hard region of pr^l GeV, believed to be properly described by 
QCD. However at that time, the exponent index obtained from fits to px 
data was equal to n~8 (or bigger), far away from the expected point inter- 
action value of n~4 |17l [T8t I19j . As one can see from our fits, this region 
of dominance of truly point-like hard interactions is still far away (albeit n 
diminishes noticeably between 900 GeV and 7 TeV). Actually, comparing 
our results with compilation of results at lower energies provided by |20] 
one observes that the naive counting rule result of n = 4 — 6 [19] seems to 
be out of reach. It means therefore, that even at the highest energies, we 
do not deal with the point-like objects expected from the naive field theory 
expectations and there is always an additional (to that resulting from a hard 
collision) pT transfer, perhaps preceded by a kind of the multiple scattering 
process or constituent scattering which make the finally observed spectra 
softer than naively expected. 

It is of interest to note in this connection that the power of n ~ 6.6 — 7.6 
in the transverse momentum spectrum at these high energies may be re- 
lated to the constituent interchange model of Blankenbecler and Brodsky 
and Gunion [17J. In the basic quark models diagrams, the power index in px 
dependence of the inclusive spectra can be inferred from the counting rule in- 
volving the collision of the active constituents [171118] (for a review, see |14]). 
If one assumes that the dominant basic high-pT process in pp — )• vrX comes 
from qq — t- qq, then the counting rule gives a transverse momentum depen- 
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dence of with n = 4, which differ from the observed power index, 

as we mentioned earher. On the other hand, if one assumes that the basic 
process is g+meson— T-g+meson, then the counting rule gives n = 8 which is 
close to the observed value. The basic process of g+meson— T-q'+meson may 
appear dominant because of the strong quark-hadron coupling. 

We close with the following observation. From what was shown here it 
seems that phenomenologically the two-parameter, QCD-inspired formula 
([5]) is as good as the two-parameter Tsallis formula with n — )• l/{q — l). The 
only difference is in the interpretation, i.e., in the possible thermodynamical 
origin (among others) of the Tsallis formula [8l[9l[10]. In this case Tsallis fits 
would cover the whole energy range of experiments uniformly interpreted 
in terms of thermal (extensive or nonextensive) model. That this view is 
reasonable was shown in papers explicitly demonstrating that nonextensive- 
thermodynamics satisfies all demands of the usual thermodynamics applied 
to systems that posses intrinsic fluctuations, memory effects, are limited 
and/or nonhomogeneous etc. [211 [22l [23] . Whether the recent CMS result 
presented here fits to this picture or rather calls for some novel explanation 
of Tsallis formula (and parameter q) remains, however, for a time being an 
open question. 
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